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Abstract 



A rigorous extraction of the deuteron charge form factors from tensor polarization 
data in elastic electron-deuteron scattering, at given values of the 4-momentum trans- 
fer, is presented. Then the world data for elastic electron-deuteron scattering is used 
to parameterize, in three different ways, the three electromagnetic form factors of the 
deuteron in the 4-momentum transfer range 0-7 fm _1 . This procedure is made pos- 
sible with the advent of recent polarization measurements. The parameterizations 
allow a phenomenological characterization of the deuteron electromagnetic struc- 
ture. They can be used to remove ambiguities in the form factors extraction from 
future polarization data. 

PACS numbers: 21.45,+v, 25.30.Bf, 27.10.+h, 13.40.Gp 
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1 Introduction 

The deuteron, as the only two-nucleon bound state, has been the subject of many theoretical 
and experimental investigations. Since it has spin 1, its electromagnetic structure is described 
by three form factors, charge monopole Gc, charge quadrupole Gq and magnetic dipole Gm, 
assuming P- and T-invariance. Measurements of elastic electron deuteron scattering observables 
provide quadratic combinations of these form factors. Since most of the data available come from 
differential cross section measurements, it has been customary, both in the data presentation and 
in the comparison with theoretical models, to use the two structure functions A and B defined 
hereafter, extracted from the cross section data by a Rosenbluth separation [jl]]. With the advent 
of tensor polarimeters and tensor polarized internal targets, polarization observables have been 
measured as well, which allow the separation of the two charge form factors. 

The purpose of this work is twofold. First, in Sect. 2, the calculation of Gc and Gq, at given 
values of the 4-momentum transfer Q, from polarization data together with (interpolated) A and 
B data is reexamined and updated with respect to previous work. 

Then, in Sect. 3, parameterizations of the three deuteron form factors, in the 4-momentum 
transfer range Q = — 7 fm _1 , are provided. Above 7 fm _1 , only small angle cross section data 
are available, preventing the separate determination of the three form factors. We have determined 
the three deuteron electromagnetic form factors by fitting directly the measured differential cross 
section [p| — p0|] and polarization []2T|-|2^] observables. This procedure eliminates the need for an 



intermediate determination of A and B, and results in a more realistic evaluation of errors for the 
form factors. 

One parameterization is used for a determination of the node of the charge form factor Gc, 



while the application of the work of Ref. [ p0| ] allows the determination of reduced form factors in a 
helicity basis. The accuracy in the determination of these form factors is limited by the assumption 
of a one-photon exchange mechanism in the first order Born approximation at low Q, and by the 
accuracy of the data at intermediate to high momentum transfers. A third parameterization was 
recently applied for a precise determination of the rms-charge radius of the deuteron [TJTJ]. At 
low Q, Coulomb distortion was taken into account to extract precise values of Gc- Applying 



this correction resolved an old discrepancy between the deuteron radius determined via (e, e') and 
N-N scattering [E2p. In the intermediate to high Q -range, other corrections such as the double 



scattering contribution to two photon exchange Q33| ] should be considered, but they are at present 
neither accurately calculated nor experimentally determined. 

2 Observables and form factors 
2.1 e-d observables 

Assuming single photon exchange, the electron-deuteron unpolarized elastic differential cross 
section can be written as 



da 

&NS " 



GUQ 2 ) + lv 2 Gl(Q 2 ) + lr,s-\Q\e e )GlM 



a NS ■ S, (1) 



dtt 

where a NS is the Mott differential cross section multiplied by the deuteron recoil factor, 9 e the 
electron scattering angle, rj = Q 2 /4M%, M d the deuteron mass; e = [1 + 2(1 + 77) tan 2 (0 e /2)] _1 
is related to the virtual photon polarization. The quantity S = A + B tan 2 (0 e /2) defines the usual 
A and B elastic structure functions. 

The tensor polarization observables t 2q , or equivalently the analyzing powers T 2q , have been 
measured as well. Their expression as a function of the three form factors, still in the one-photon 
exchange approximation, is given by: 

- V2 ■ S ■ ho = \vG c G Q + % - V 2 Gl + Ive-'Gl (2) 

V3 ■ S ■ t 21 = 2 V L + V 2 sin 2 °A ' G m Gq sec | (3) 
-2V3-S-t 22 = V G 2 M . (4) 
2.2 Calculation of G c and Gq 

The charge form factors are here extracted from t 2 o(Q,9 e ) data, together with A(Q) and 
B(Q) (interpolated) data. The analyses presented in [ |2~2"| , [2~6T ] need to be updated, because of new 
^20 [EU-EH-EIB anc ^ ^ ©Si ^ ata - m particular, the parameterization of A used in [ ^6] ] gave a very 
small weight to the then only existing high Q data [|5|] and is lower than the new data [QQ] around 
4.5 fm _1 . Furthermore, we present here a more compact solution and a more rigorous treatment 
of errors. 



For our purpose, it is useful to define new quantities A = A — B/2(l + rj) and £20 Q2q], derived 
respectively from A and t 2 o by eliminating the magnetic contribution: 

\ V G C G Q + l V 2 G 2 Q _S-t 20 + B/4y/2e(l + 77) 



V2 {<% + ir?(? Q ) " A 
Using the reduced form factors g c = G c /\fA) and g Q = 2t]Gq/3-\/A^, lead to: 



(5) 



0c + 2$ = l (6) 
2^q + 9 Q = V = -ho/ V2 (7) 

where p (or conventionnally pzz) is the tensor polarization in Cartesian notation (also called align- 
ment). There are four solutions to these equations given by 

with A = 8(1 — p){\ + p) and from ([7]). The physical solution is easily selected at small Q 
from the static moments (<7c(0) = 1> 9q(Q) — 0). It corresponds to the choice of a minus sign 
in (|8|) and of qq > 0. Since t 2 o an d hi, both proportional to Gq, do not cross zero at a same 
value of Q $n$m , gQ has to remain positive over the whole range considered in this work. The 
two remaining solutions (gn 7 gc) and (gn,ga) cross each other at values Q m in and Q m ax where 
i 2 Q reaches its extrema —y/2 and +\/2/2 (A = 0). The physical solution must switch from "— " 
to "+" at Q = Qmin and then back to "— " at Q ma x in order to ensure a continuity of the form 
factor derivatives. For polarization data close to these extrema, Q may be below or above the a 
priori unknown Q m in (or Q ma x), and the choice of solution is ambiguous. Qmin, from our three 
global fits to the e — d data (see Sect. 3), is determined to be close to 3.3 fm -1 . On the other hand, 
there are not enough polarization data to constrain the value of Qmax, so that the above mentioned 
ambiguity remains around Q ~ 6 — 8 fm -1 . This is the case for the two points at highest Q in 

An additional complication arises for five polarization data points in Refs. [ffH- P3|JZ5| , p7| ] which 



lay partially outside the physical region — y2 < £20 < 1/ V%- This situation is quite probable for 



points with finite errors close to a physical limit [|37|]. For the sake of extracting Gc and Gq, the 



interval of 68.3% confidence level [i 20 — At 20 , ho + Ai 20 ], and eventually the most probable value 
t 2 o, are then modified according to the method presented in [|38|]. The resulting confidence interval 
is entirely within the physical region (A > 0). In this particular case, the modified values of p are 
used in (QJ^) instead of the measured ones. As a result of this procedure, the errors on the form 
factors may be asymmetric. 

The calculated values of Gc and Gq, corresponding to all measurements of ho, are presented 
in Table 1 and Fig. 1. The later also shows results of parameterizations to be discussed in Sect. 
3. Uncertainties come from the quoted errors in ho, combined quadratically with errors on A and 
B reflecting the spread of the data (for example, at 5 fm _1 , 8.5 and 17 % respectively). For the 
two points of highest Q, the two solutions of are given. The first one is preferred, based 
on theoretical guidance and on the parameterizations discussed below. Only parameterization I 
(Sect. 3.1) favors the second solution for the point at Q = 6.64 fm _1 . Note that t 20 need not 
necessarily reach its maximum allowed value, in which case the first ("+") solution would prevail 
from Q = Q min up to the undetermined node of Gq, or to the second minimum of t 20 , whichever 
occurs first. 

3 Parameterization of the form factors 

The three paramaterizations described below are determined through a x 2 minimization in- 
volving 269 cross section data points P-|20|] and 39 polarization data points [|2T]— 129|] . In most 



polarization data, and in some cross section data, the systematic uncertainties are dominant and 
may vary from point to point in a given experiment. The error considered in the x 2 minimization 
is then the quadratic sum of the statistical and systematic uncertainties. The uncertainties on the 
parameters are given by the error matrix. For data where an overall normalization uncertainty 
may apply, the resulting systematic uncertainty of the fitted parameters have been evaluated by 
changing each individual data set by the quoted error and re-fitting the complete data set. This last 
procedure was carried on only with parameterization III (Sect. 3.3). 

The x 2 P er degree of freedom {x 2 l^d.f) all exceed the value of 1, because of systematic 
differences between some data sets, at the limit or beyond the quoted systematic uncertainties. 
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Among the most recent experiments, this is the case for the A measurements of Refs. [g^J], and in a 
lesser extent for the t 2 o measurements of Refs. [ [2~T| , [2"6T ] . The fits then give an average representation 
of the data, though biased toward experiments with a larger number of data points. 
3.1 Parameterization I 

In the first parameterization (I), each form factor is given by: 



G X (Q 2 ) = G x (0) 



1 



Q_ 



X . 



1 + ^2a Xi Q 



2 1 



i=l 



(9) 



with X = C, Q or M. This expression has the advantage of displaying explicitly the first node 
Q x of each form factor. The normalizing factors Gx(0) are fixed by the deuteron static moments. 
With 18 free parameters, a fit is obtained with x 2 /Nd.f. = 1.5. 
3.2 Parameterization II 

Another parameterization (II) has been proposed by Kobushkin and Syamtomov [00]. Each 
form factor is proportional to the square of a dipole nucleon form factor G d and to a linear com- 
bination of reduced helicity transition amplitudes go,gi,g2- 



1 ' r ^ 

Gq 
\ Gm j 



G 



M{q) 



( \ 

go 

9i 
K 92j 



(10) 



Each of these amplitudes is parameterized as a sum of four Lorentzian factors: 

i 



9k = Q k Y.' 



a ki 



(ID 



For each k, the d\ i follow an arithmetical suite defined by 2 independent parameters. In addition, 
an asymptotic behavior dictated by quark counting rules and helicity rules valid in perturbative 
quantum chromodynamics (pQCD), together with the normalization conditions at Q = 0, imply 
6 relations between the parameters and p^]. As a result, each amplitude is described by 
4 independent parameters. New parameters are obtained here, due on one hand to a newer data 
base, and on the other hand to the fitting of the differential cross sections instead of A and B. 
With 12 free parameters, a fit to the data set is obtained with x 2 /Nd.f. = 1.8, whereas the original 
values of the parameters in Ref. [ |30| ] yield x 2 /Nd.f. = 7.5. This parameterization, in contrast with 



the two other ones presented in this paper, can be extrapolated well above 7 fm -1 , albeit with 
some theoretical prejudice. We confirm the observation of Refs. [ J5U| , P4| that the double helicity 
flip transition amplitude g 2 has a magnitude comparable to the zero helicity flip amplitude g in 
the Q-range considered here, which means that these amplitudes are not in the asymptotic regime 
expected from pQCD. 

3.3 Parameterization III 



The third parameterization (III) employs a Sum-of-Gaussians (SOG) Q35Q. The form factors 
are written as 

Although our interest here lies in its Q-space version, the parameterization is better described in 
configuration space where it corresponds to a density p(R) written as a sum of Gaussians placed at 
arbitrary radii Ri, with amplitudes A4 fitted to the data, and a fixed width 7. The distance R refers 
to the distance of the nucleons to the deuteron center of mass. The parameterization represents a 
totally general basis and the following applied restrictions are justified on physics grounds. First, 
one does not expect structures smaller than the size of the nucleon, which determines the width 



7 to be the size of the proton (7^3/2 = 0.8 fm). Second, the spacing between Gaussians is 
chosen slightly smaller than this width: 0.4 fm or 0.5 fm. Third, the Gaussians are placed at 
radii R L < R max = 10 fm, which is justified given the fact that one can easily specify the radius 
at which the tails of densities give no significant (< 10 -3 ) contribution to Gx(Q)- In addition, 
outside the range of the NN-force, the deuteron wave functions have an analytic form which is 
well known and depends only on the deuteron binding energy. Thus, for radii Ri > 4 fm, one 
can impose this shape and fix the ratio of the amplitudes A*. Each form factor is then determined 
with 1 1 free parameters: 10 Gaussian amplitudes A 1 to A 10 , corresponding to Ri < 4 fm, and one 
overall amplitude for the shape-given tail at R > 4 fm. With a total of 33 independent parameters, 
a x 2 /Nd.f. of 1.5 is obtained in the fit. 
3.4 Results and discussion 

The resulting form factors from the three parameterizations are shown in Fig. 1. As functions 
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of two variables (Q and 6 e ), the fitted quantities cannot be easily represented together with the pa- 
rameterizations. In order to illustrate the quality of the fits, we present plots of relative differences 
of A and B, and of £20 (Q) in Fig. 2. t 2 i and t 22 are equally well fitted, which constitutes, within 
experimental uncertainties, an indication of the coherence of equations ([T^^l]), and therefore of 
the consistency of the one-photon exchange approximation. 

From the average and dispersion between the three parameterizations, combined with the fit 
uncertainty on Qq, the node of the charge form factor is determined to be located at 4.21 ± 0.08 



fin -1 , a value governed by the t 20 results of Refs. [ j2l"1 , [2~6| ] . Assuming as we do here implicitly that 
these two data sets have the same weight, the location of this node is not quite consistent with a 
relation between the two- and three-nucleon isoscalar charge form factors, established with various 
N — N potentials Q36[l. The secondary maximum of \Gq\ is very flat, so that its location (5.3 ± .5 
fm -1 ) is not determined very precisely. Its magnitude (.0038 ± .0003) is clearly inconsistent with 
the corresponding one of the three-nucleon isoscalar charge form factor, still within the same 
model calculations The £ 21 results of Ref. pT|], though of limited accuracy, help confirm a 



node of the magnetic form factor \^ at 7.2 ± 0.3 fm -1 . As for the first node of Gq, according 
to most theoretical models, it should appear at a higher value of Q, above the range where our 
parameterization method applies. The value Qq = 7.7 ± 0.6 fm -1 given by parameterization I is 
probably the smallest possible value allowed by the present data. It is due to this parameterization 
following the downward trend of the £ 20 data point at the highest Q (see Fig. 2). This trend 
however is not statistically significant. Parameterization II, when extrapolated, suggests a much 
higher value of Q for the node of Gq. Finally, from 



2 a 

r = — 6 



= 6 

Q 2 =0 



aci + (Q° c )- 2 , (13) 



dQ 2 

we calculate the root mean square charge radius of the deuteron r = 2.094 ± 0.003 (stat.) ±0.009 
(syst.) fm. The statistical uncertainty is given by the error matrix from parameterization I, while 
the systematic uncertainty is evaluated with parameterization III (see above remark about normal- 
ization uncertainties on individual data sets). This radius is 1.7% smaller than the value r = 2.130 



fm reported in Q31|], consistent with expectations in the absence of corrections due to Coulomb 



distortion. 



4 Conclusion 

The extraction of the charge form factors G c and Gq from experiment, at given values of 
Q, has been reexamined. The solutions were expressed in the most compact and physical way, 
while a new treatment of errors was applied to polarization data at or beyond physical limits. 
The existing electron-deuteron elastic scattering data were used for direct parameterizations of 
the three deuteron electromagnetic form factors, up to Q = 7 fm -1 . The numerical results may 
be requested from the authors^ and will be updated as new data become available in the future. 
The inferred value of Q m in — 3.3 fm -1 corresponding to the minimum of t 2 o could be used, or 



recalculated with such global fits, for future experiments in this Q-range [ p9l|40[ ], in order to resolve 
the discussed ambiguities in the form factors calculation. These future experiments should help 
confirm, or adjust, the exact value of the node of the charge form factor: this location is sensitive to 
the strength of the N — N repulsive core, to the size of the isoscalar meson exchange contributions 
and to relativistic corrections. The observation of the node of the magnetic form factor [|8|j21|] 
should be confirmed in a more precise experiment pT]]. Together with the determination of the 
secondary maximum of \Gc\ [p]], this would complete the full characterization of the deuteron 
electromagnetic structure up to Q — 7 fm -1 . 
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TABLES 



Q 

(fm' 1 ) 


%(70°) 


^20 


G c 


G Q 


Ref. 


0.86 


-.30 (±.16) 


-.30 (±.16) 


.627 (±.011) 


47. (±25.) 


[@] 


1.15 


-.181 (±.070) 


-.178 (±.071) 


.474 (±.008) 


12.0 (±4.7) 


Hi 


1.58 


-.400 (±.037) 


-.402 (±.038) 


.289 (±.006) 


8.66 (±.81) 


d] 


1.74 


-.420 (±.060) 


-.423 (±.063) 


.238 (±.005) 


6.19 (±.90) 


ill 


2.026 


-.713 (±.090) 


-.734 (±.095) 


.160 (±.005) 


5.51 (±.73) 


[@i 


2.03 


-.590 (±.130) 


-.604 (±.138) 


.163 (±.005) 


4.50 (±1.02) 


i& 


2.352 


-.896 (±.093) 


-.945 (±.101) 


.100 (±.004) 


3.49 (±.41) 


[0 


2.49 


-.751 (±.153) 


-.792 (±.169) 


.087 (±.004) 


2.17 (±.48) 


[0] 


2.788 


-1.334 (±.233) 


-1.473 (±.267) 


3.71(lJ;g)xlO- 2 


2.59 (±.073) 




2.93 


-1.255 (±.299) 


-1.401 (±.347) 


3.45 (™)xl0- 2 


1-85 (t;ll) 




3.566 


-1.87 (±1.04) 


-2.20 (±1.26) 


1.53 (™)xl0- 2 


.651 (+:&S) 




3.78 


-1.278 (±.186) 


-1.476 (±.228) 


1.25 (t^)xlO- 2 


.474 ct:S3) 


[§] 


4.09 


-.534 (±.163) 


-.567 (±.193) 


-1.14 (±1.6)xl0~ 3 


.383 (±.015) 


[0] 


4.22 


-.833 (±.153) 


-.913 (±.179) 


1.63 (t^)xlO- 3 


.325 (±.013) 


[§] 


4.46 


-.324 (±.089) 


-.320 (±.100) 


-2.39 (±.61)xl0- 3 


.245 (±.010) 




4.62 


-.411 (±.187) 


-.417 (±.207) 


-1.63 (±1.14)xl0- 3 


.208 (±.009) 




5.09 


.178 (±.053) 


.208 ±(.056) 


-3.87 (±0.30)xl0- 3 


.119 (±.006) 




5.47 


.292 (±.073) 


.312 (±.075) 


-3.48 (±0.32) xl0~ 3 


.080 (±.004) 




6.15 


.621 (±.168) 


.630 (±.170) 


-3.19 (±0.55)xl0~ 3 
-4.20 (t-j 2 )xl0- 3 


.034 (i;°°l) 
.019 (±.007) 




6.64 


.476 (±.189) 


.478 (±.189) 


-1.89 (±0.38)xl0~ 3 
-3.13(i; 2 ^)xl0- 3 


.023 (±;™ 2 ) 
.008 (±.004) 
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TABLE I. Calculated values of t2o(70°), £20, Gc and Gq corresponding to all £20 measurements. In 
parantheses, statistical and systematic uncertainties are added in quadrature. For the last two points, the two 
solutions are given (see text). 
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FIGURES 



FIG. 1. 

Deuteron form factors Gq, Gq and Gm as a function of Q. The data for Gc and Gq are from 
Table 1, corresponding to t 2 o measurements of Refs jnj ] (solid diamonds, and open diamonds for 
the second solution), [£2p (star), [fZ3]] (open squares), P^j29| ] (triangles up), p5| ] (open circle), 



[ P6| ] (full squares), [ |27| ] (triangles down), [[28]] (full circles). The Gm data corresponds to the B 
measurements of Refs. [§]] (open diamonds), [|8|] (open circles), [ |T0| ] (stars), [ |20l ] (full circles). 
The curves are from our parameterizations I (solid line), II (dot-dashed line) and III (short dashed 
line). 



FIG. 2. 

(a) AA/A, in %: deviation of A with respect to parameterization I, arbitrarily taken as a 
reference line; for clarity only the data from Refs [g] (full diamonds), [Q] (full circles), [|5p (open 
circles), [ ]T^[ ] (triangles), [|TS]j (open diamonds) are reported, (b) AB/B, in %. (c) t 2 o, with 
physical domain delimited by dotted lines. For B and t 2 o data legend, as well as curves legend, 
see Fig. 1. 
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